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ABSTRACT 

We consider the radial migration of vortices in two-dimensional isothermal gaseous disks. We find 
that a vortex core, orbiting at the local gas velocity, induces velocity perturbations that propagate 
away from the vortex as density waves. The resulting spiral wave pattern is reminiscent of an embedded 
planet. There are two main causes for asymmetries in these wakes: geometrical effects tend to favor 
the outer wave, while a radial vortensity gradient leads to an asymmetric vortex core, which favors the 
wave at the side that has the lowest density. In the case of asymmetric waves, which we always find 
except for a disk of constant pressure, there is a net exchange of angular momentum between the vortex 
and the surrounding disk, which leads to orbital migration of the vortex. Numerical hydrodynamical 
simulations show that this migration can be very rapid, on a time scale of a few thousand orbits, 
for vortices with a size comparable to the scale height of the disk. We discuss the possible effects of 
vortex migration on planet formation scenarios. 

Subject headings: accretion disks — hydrodynamics — planets and satellites: formation — waves 



1. INTRODUCTION 

Planets are thought to form in circumstellar disks of 
gas and dust that are commonly found around young 
stars. Within these disks, planet formation requires that 
solid material has to grow 13 orders of magnitude in size, 
from sub-micron sized interstellar grains all the way to- 
wards terrestrial planets. Large-scale, long-lived vortices 
can play an important role in this process, since they are 
very efficient in collectin g solid particles in their center 
( Barge fe Sommeria|1995[ ). It has been proposed that this 
can speed up the formation of protoplanets enormously 



(Lyra et al. 2009). It has also been claimed that vor- 



tices can e ffectivel y transport a ngular momentum out- 



ward ( Johnson & Gammic 2005 ) , allowing for an inward 
accretion now m disks that have too low ionization frac- 
tions to be unstable to the MRI (magneto- rotational in- 
stability, |Balbus fc Hawley||1991| ). 

The emergence and subsequent survival of vortices in 
protoplanetary disks has received a lot of attention re- 
centl y. They can appear due to Rossby wave instabili- 
ties ( Lovelace et al.|1999l, edge instabilities of planetary 
gapsl Kollcr et al. 2003 foe Val-Borro et al.]|2006l [20071 
Lin fc Papaloizou 2010), in 3D circulation mo dels ([Bar ' 
;us|2[)J5 ) and in MHD turbulence ( Fromang 



ranco & Marcusi: 



fc INelson 2005 JT Vortex generation due to a baroclinic 
instability in dis ks with a radial entropy gradi ent was 



first proposed in Klahr & Bodenheimer (2003). How- 
ever, this study was hampered by numerical issues, and 
coul d subsequently not be repro d uced in a local shearing 
box ( Johnson fc Gammie||2006 ) . Petersen et al. (2007a I 
pointed out the importance of thermal diffusion in this 
problem, and showed, using anelastic global si mulations 



that vortices can be produced and maintained (Petersen 
et al.|200"7b ) in disks with a negative radial entr opy gra- 



dient. More recently, Lesur & Papaloizou (2010) showed 
that the baroclinic instability is in tact subcritical and 
needs a finite amplitude disturbance to get going. They 
also showed that the vortices produced from this Sub- 
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critical Baroclinic Instability (SBI) can survive in 3D, 
despite the fact that they are unstab le to parametric in- 
stabilities ( Lesur fc Papaloizou||2009[ ) . 

In this paper, we take the emergence of vortices to be 
given, and focus on one aspect of their subsequent evolu- 
tion that has not been studied so far: orbital migration. 
It is well known that solid bodies embedded in Keplerian 
gas disks change their orbit through interaction with the 
gas. Small bodies experience a head wind from the gas, 
which is partly pressure supported, caus ing them to lose 
angul ar momentum and spiral inward ( Weidenschilling 
1977[ ). The gravitational perturbation of a planet em- 
bedded in t he disk leads to wave excitat ion at Lindblad 
resonances dGoldreich &JlY emaine|1979| , which is asym- 
metric in general ([Ward 1997j|, lead ing to inward migra - 
tion for reasonable disk parameters ( Tanaka et al.| 2002). 
This is called Type I migration. Interestingly, vortices 



lations (see Lesur & Papaloizou 


2010 


Heinemann & Pa- 


paloizou 


2009a|b|), and one mig 


tit expect that the same 



etary Type I migration could lead to radial migration of 
vortices. This is the subject of this paper. 

The plan of the paper is as follows. We introduce the 
disk models in Sect. [2j and discuss the numerical meth- 
ods used in Sect. [3j We then look at a particular example 
of vortex migration in Sect. HI showing that this migra- 
tion is due to the emission of density waves into the disk. 
We discuss the various physical mechanisms that lead 
to asymmetric wave emission in Sect. [5] show that this 
leads to vortex migration in Sect. [6j and compare these 
findings to numerical results in Sect. [7] We discuss our 
results in Sect. [8] 

2. BASIC EQUATIONS AND DISK MODELS 

We consider thin disks in the two-dimensional approx- 
imation, and work with vertically averaged quantities 
only. The governing equations are then the continuity 
equation 

— -i- v . t,v = n m 
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where £ is the surface density and v = (v r ,r£l) T is the 
two-dimensional velocity in cylindrical coordinates (here 
T indicates transpose), and the momentum equation 



v • Vv 



(2) 



with p the vertically integrated pressure and $ the grav- 
itational potential. We consider inviscid disks only, and 
ignore the self-gravity of the gas, so that $ is due to the 
central star only. We will consider a strictly isothermal 
equation of state, p — c^S, with c s the sound speed. 

The disk surface density profile is a power law initially, 
X = £ (r/r ) -Q , with r a reference radius, most often 
the initial location of the vortex. Since we do not con- 
sider self-gravity, So is arbitrary. The sound speed is 
set by choosing a scale height H at ro, Hq, which de- 
termines the sound speed c s = Hq^Iq. We then have 
H = Ho(r/ro) 3 / 2 . The angular velocity SI is Keplerian, 
initially, with a small correction due to a radial pressure 
gradient. In some of the analysis, it is advantageous to 
have a constant aspect ratio H/r. This can be achieved 
by having a radially varying, but constant in time, sound 
speed. Physically, this corresponds to a disk that relaxes 
to an equilibrium temperature profile on a very short 
time scale (the cooling time scale is formally zero, so no 
temperature fluctuations are allowed). 

It is possible to combine equations ([!]) and Q into a 

single equation for the vortensity w/S, where uj = k • 
(V x v): 



d_ 

Ft 



(I) 



v V 



VS x Vp 

E3 ' 



(3) 



showing that in 2D barotropic flow, for which p ~ p(S), 
vortensity is conserved along streamlines. This of course 
includes the isothermal case. It is no longer true when 
the sound speed varies radially, which makes the flow 
non-barotropic, or when 3D motions are allowed for. In 
that case, a more general quantity (V x v) • VQ/p, with 
Q any quantity that is conserved on fluid elements and p 
the density, is conserved along streamlines. For example, 
for an adiabatic flow Q can be taken to be the specific 
entropy. The resulting quantity is sometimes called the 
potential vorticity. 

On top of the equilibrium configuration we introduce 
a vortical perturbation. In global models, this is done by 
applying a circular velocity perturbation with a specified 
maximum over a circular patch of the disk. The velocity 
perturbation is usually a sizeable fraction (~ 0.5) of the 
local sound speed, and the patch a sizeable fraction of the 
local scale height, typically Ho/2. The velocity pertur- 
bation tapers off exponentially away from the maximum. 
The exact form of the perturbation was found not to be 
important, as the system relaxes to a similar vortensity 
profile on a dynamical time scale regardless. Note that 
the subcritical baroclinic instability will tend to produce 
vortic es of a size comparable to H (Lesur & Papaloizou 
2010|. 



3. NUMERICAL METHODS 



For our global disk models, we solve equations ([I]) 
and ([2]) in cylindrical geomet ry (r, ip) using rodeo 
( Paardekooper fc Mellema|2006 ), which is a second-order 
finite volume method using an approximate Riemann 
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Fig. 1. — Relative perturbation of vortensity (top) and density 
(bottom) after 10 orbits at r = 1 for an isothermal disk with Ho = 
O.lrrj, a = 3/2, and an initial velocity perturbation of 0.5c s over a 
circular region of radius Hg/2 around r = X,tp = 7T. 



solver due to Roe (19811. It has been success fully ap- 



plied to the study of disk-embedded planets (|de Val-| 
Borro et al.|2006 Paardekoop er &; Papaloizou 200 81 and 
disks in tight binary systems (|Paardckoo per et al.|2 008 ) . 

The computational domain consists of an annulus at 
a reference radius ro, with inner and outer boundaries 
located at 0.2 r and 2.5 r , re spectively. W e use non- 
reflective boundary conditions ( Godon| |l996), to allow 
the waves excited by the vortex to leave the compu- 
tational domain freely. Although designed to be non- 
reflective for ID simple waves only, we have observed 
no wave reflection in our 2D non-linear simulations. As 
explai ned in more detail in |Paa rdckoop er fc Mellema| 
( 2006 ) , the use of a characteristics-based Riemann solver 
ensures that information is always drawn from the right 
places. At the boundary, any incoming characteristics 
are treated as if they originate from an unperturbed 
disk. This ensures that no waves enter the computa- 
tional domain. We consider the full 2n in azimuth. For 
Hq = O.lro, we use a typical grid size of 1536 in the 
radial and 6144 in the azimuthal direction, varying this 
up and down by factors of 2 to study effects of resolu- 
tion. Our typical grid therefore has 67 zones per scale 
height. When considering smaller values of Hq, we in- 
crease the resolution to maintain a fixed number of grid 
cells per scale height. Within the code, we set ro = 1 for 
numerical convenience. 

4. A FIDUCIAL CASE 
4.1. Orbital migration 



3 



In this section, we consider an isothermal disk with 
a = 3/2 and Hq = O.lro. Note that this is a disk with 
constant vortensity, and that, for numerical convenience, 
it is relatively thick for a protoplanetary disk. The de- 
pendence on H is discussed in Sect. |7j We give the 
disk a circular velocity perturbation of magnitude c s /2 
over a scale Ho/2. The resulting vortensity and density 
perturbations after 10 orbits at r — 1 are shown in Fig. 

From Fig. [T]we see that the vortensity perturbation re- 
mains localized, and that its radial extent roughly corre- 
sponds to the size of the initial perturbation. The density 
perturbation inside the vortex is modest (~ 5%) com- 
pared to the vortensity perturbation (~ 300%), which in- 
dicates that on scales <C H , the disk response to the per- 
turbation is approximately incompressible, as expected. 
For a vortex in equilibrium in a shear flow, there exists 
a rela tion between the vorticity perturbation and aspect 
ratio ( jKida|ll981[ |Lesur k. Papaloizou||2009| . In equilib- 
rium, the relative perturbation in vorticity m a Keplcrian 
disk should be {lo — uj )/uj n = — 3(x+ l)/(x 2 — x)i where 
X is the aspect ratio of the vortex. The strong vortex 
depicted in Fig. [T] has an aspect ratio of ~ 2.5, which 
would lead to an equilibrium relative vorticity perturba- 
tion of ~ —3. Noting that we can neglect density pertur- 
bations in the almost incompressible disk response, and 
can therefore equate vorticity and vortensity relative per- 
turbations, this is in good agreement with the top panel 
of Fig. [T] We note that this equilibrium is established 
for a range of vortex aspect ratios within a few orbits 
of the vortex. We therefore conclude that the vortices 
quickly adopt an equilibrium configuration, regardless of 
the initial perturbation. 

The most notable density features start at a radial dis- 
tance of approximately H from the vortex, where the 
flow relative to the vortex becomes supersonic and den- 
sity waves can be excited. The resulting two-armed spiral 
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Fig. 2. — Time evolution of the radial location of the vortex, 
for three different resolutions. Disc and vortex parameters are the 
same as in Fig. [l] 



gration is an important process to consider. Additional 
runs that included viscosity, with the kinematic viscos- 
ity v parametrised using the a-prescription, v — a v c s H, 
indicated that the numerical dissipation of these vortices 
roughly correspond to a v — 10~ 4 , 10 -5 and 10~ 6 re- 
spectively for the low, medium and high resolution cases 
shown in Fig. [2] 

4.2. Wave action 

Density waves propagate inwards and outwards away 
from the vortex. Provided the vortex is not too strong, 
we expect the waves to be in the linear regime close to 
the vortex, becoming no n linear at larger distances w here 



wave pattern is reminiscent of an embedded planet (|Bry-| shocks are formed (e.g. |Goodman fc Rafikov||2001[ ) and 



den et al. 1999), with one inner wave propagating into 
the inner disk, and one outer wave propagating into the 
outer disk. For embedded planets, it is well-known that 
geometrical effects tend to favor the outer wave, which 
removes angular m omentum from the p lanet and causes 
inward migration (Tanaka et al. 2002). An additional 
complication in the vortex case is that the vortex may 
be asymmetric with re spec t to rrj. We deal with vortex 
asymmetries in Sect. |5.4[ here, we just note that for 
a = 3/2, the vortex is symmetric to a high degree. The 
resulting migration of the vortex is shown in Fig. [2] 

Indeed the vortex migrates inward. Results are shown 
for three resolutions in Fig. [2j The lowest resolution has 
16 cells per scale height, radially, at ro- In other words, 
the vortex is resolved by approximately 8 cells only, in 
the radial direction. In this case, the vortex weakens 
through numerical diffusion, which starts to slow down 
the migration rate after 10 orbits. Doubling the reso- 
lution pushes this time towards 20 orbits, and doubling 
it again gives a steady migration rate for at least 40 or- 
bits. At early times, the migration rate is similar for all 
resolutions, and we conclude that this migration rate is 
converged with respect to numerical resolution. It results 
in a migration time scale of ro/\r\ s» 2000 fig , or 300 
orbits at r = ro. It is clear that for this disk, vortex mi- 



their amplitude reduces. While in the linear regime, the 
total rate of flow of angular momentum across a radial 
location that is associated with the waves, or equiva- 
lently the wave action, is conserved. When the non linear 
regime is entered and dissipation occurs, the wave action 
decreases. Angular momentum is advected inwards in 
the radial direction at a rate given by 



A = 2i 



(TiVrV 



r u ip I 



(4) 



where the angle brackets denote an azimuthal average. 
Because there is no net mass flow associated with linear 
waves, the azimuthal velocity perturbation alone may be 
used in Q in this case. We remark that, by making 
use of thelinearized equations governing the waves is is 
possible to write Q in a different form. The waves have 
a pattern speed fJvort which is also the angular velocity 
of the vortex. Thus for the perturbations associated with 
them we have 

d d 

Qvort ^ = at' (5) 

Accordingly the perturbed azimuthal component of the 
equation of motion gives 
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Fig. 3. — Wave action after 10 orbits for the same vortices as in 
Fig. [2] 

where E' is the surface density perturbation. Noting that 
to linear order, the Lagrangian displacement £ r satisfies 



(Q — Q VO rt) 



dtp 



we see that 

(n 



n 



vort) ( w ¥> + 2f2^ r 



rE ■ 



(7) 



(8) 



Using ^ to substitute for v v in Q and making use of 
Q , it is readily seen that we can rewrite Q in the form 



A = 2ir 



rc 2 s Ylv r 



(9) 



We can now verify that this migration is due to the 
density waves by studying the behavior of wave action. 
The change in the angular momentum content of the fluid 
on account of the presence of the vortex is 



J (Er 2 ft 



E^r Cli)dS, 



(10) 



where subscripts i denote the initial conditions, and the 
integral is taken over the whole vortex surface area. We 
expect dL vm t/dt — AA, with A A the difference in wave 
action defined above evaluated for the outer and the inner 



wave. 



We show the wave action for three resolutions in Fig. 
[3| Equation [9] only has meaning in the wave region, at 
distances of at least 2H/3 from the vortex. Although our 
lowest resolution run clearly does not resolve the waves 
well enough to get the correct wave action on both sides, 
A A is very similar to the higher resolution runs. In Sect. 
[7j we make sure our resolution is at least as high as that 
of the dashed line in Fig. [3] The similarity of the wave 
action for the highest resolution runs translates into a 
similar migration rate after 10 orbits (see Fig. [2|. We 
have measured the change in angular momentum of the 
vortex between 5 and 10 orbits, and compared this to 
what is expected from the wave action. The results are 
summarized in Table [T] 

We point out that there is a large uncertainty in the 
measured AL. This is because L vort can only be mea- 
sured to an accuracy of ~ 5% because of ambiguities in 



TABLE 1 

Vortex angular momentum and wave action 
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a Change in orbital distance of the vortex between 5 and 10 orbits. 
b Angular momentum transport by density waves. 
c Change in angular momentum of the vortex between 5 and 10 
orbits. 




Fig. 4. — Time evolution of the wave action difference, for the 
same three resolutions as in Fig. [2] 

defining the vortex as a separate entity from the back- 
ground flow, and AL results from a subtraction of two 
almost equal numbers (typically AL ps 0.1L vort ). We 
therefore estimate that there is an uncertainty of at least 
50% in the measured values of AL, which makes the 
agreement between AAAt and AL reasonable (see Ta- 
ble [I]). Therefore, the inward migration of the vortex is 
consistent with the wave torque. 

We show the time evolution of the wave action differ- 
ence, or wave torque, in Fig. |4j for the same three reso- 
lutions as in Fig. [2] In all cases, a steady torque is set 
up after 3 orbits, at which time the vortex has reached 
its equilibrium with the background shear. After the 
equilibrium has been set up, the slow time evolution of 
the torque is due to numerical diffusion (the decline in 
the low resolution case) and slowly varying background 
disk properties due to the changing radial location of the 
vortex (the slow increase in the highest resolution cases). 

5. WAVE ASYMMETRIES 

5.1. Linear equations 

In order to understand the origin of the waves emit- 
ted by accretion dis k vortices, let us int roduce the local 
shearing box model ( Hawley et al.fl995| ). In this model, 
we neglect curvature effects in equations ([l])-((2]) and in- 
troduce a local cartesian frame corotating witn the disk 
at r by defining x = r — r Q and y — r$tp. In this model, 
equations and ^ read: 



as 

dv 



+ V-Ev=0, 
+ v ■ = 



(11) 



-4 V In E - V(-gfiV) -2ilXv, 



(12) 

where 57 = O(r ) is the Keplerian frequency at r , 
q = 3/2 for a Keplerian rotation profile and we have 
assumed an isothermal gas. In the following, we consider 
a vortex located at r = tq and perturbing the surround- 
ing flow. The vortex core is a non-linear solution to the 



■5 



above equations (see e.g. Lesur fc Papaloizou|2009 ), and 
we consider the linear perturbations produced by the 
vortex core at long distance. Since the vortex is, in first 
approximation, a quasi-steady structure, we consider sta- 
tionary perturbations of the Keplerian flow v Q = —qtte y 
introducing 



T,'(x,y)=T,(x,y) - S , 
u(x,y)=v -v . 



(13) 
(14) 



As a further simplification, we only consider the case of 
a constant background density profile So- To reduce the 
problem to a system of ordinary differential equations, we 
introduce a Fourier decomposition in the y direction: 



S' = S'(a;) exp(i/cy), 
u = u(x) exp(iky). 



Using this decomposition in (11) and (12) yields 



ier(x)S' + Sq ( -j^- + ikuy 



—ia(x)u x — 2ilu v H — f-— — =0 
So dx 

S' 

—ia(x)uy + (2 — q)D,u y + ic s k— = 

S 



(15) 
(16) 



(17) 
(18) 
(19) 



where we have defined o~{x) = qftx. After some algebra, 
one gets the following equation for the azimuthal velocity 
u y : 



dx 2 



+ A 2 (x)u y = 0, 



where 



A 2 (x) = 



a 2 (x)- 2(2 -q)fl 2 



-k 2 



(20) 



(21) 



This equation describes the propagation of a linear 
acoustic-inertial wave in the shearing box. As stated 
above, we consider a nonlinear vortex core located at 
x = emitting a wave at large distance. Assuming the 
radial extent of the vortex is very small compared to the 
disk thickness H — c s /ft, although this need not be the 
case for the azimuthal extent, the vortex core can be 
seen in first approximation as a boundary condition at 
x = for the wave emission problem. Without loss of 
generality, we will consider the case of a wave emitted in 
the x > region, the x < solutions being deduced by 
symmetry arguments. 

It is worth noting that ( p0| is in fact the vorticity con- 
servation equation in a compressible medium, with a con- 
stant background vorticity (2—q)H,. In a global disk how- 
ever, background vorticity is not constant and density 
might also depend on radius. In this context, one would 
write a conservation equati on for the vortensity V X v/T, 
which would be similar to ( 20 ) if the background vorten- 
sity was constant. Therefore, the wave emission model 
we are discussing here should be understood as a local 
model for a constant vortensity disk. 

5.2. Wave solutions 

The above equation admit solutions in the form of 
parabolic cylinder functions. However, to understand the 
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Fig. 5. — Supersonic solution wave fronts as obtained from ( |2."i[ i 
assuming x s = 1. The sonic lines are indicated by dashed lines. 

physical origin of vortex wa ves, we present here WKB so- 
lutions. We first note that (20) has a turning point at 



qQk 



(22) 



where we have introduced uj a i — ^/k 2 c 2 + q(2-q)n 2 , 
the natural frequency of acoustic-inertial waves. This 
particular location corresponds to the sonic line of the 
vortex. This is the line along which the flow speed is 
equal to the wave speed as given by the phase velocity. 
In the limit of large k this phase velocity becomes equal 
to the sound speed hence the term sonic line. For x < x s 
("subsonic" region), the wave has an exponential shape 
and it corresponds to the logarithmic tail of an incom- 
pressible vortex. Using first order WKB theory in the 
limit x <C x Sl we get: 



Aexp 



A 2 (u) 



1/2 



du 



B exp 



A 2 (u) 



1/2 



du 



(23) 



On the other hand, for x > x s ("supersonic" region), the 
wave is essentially sinusoidal and can travel for very large 
distances. In the limit x 3> x s , we obtain: 




A 2 (u) 



1/2 



du 



1/2 



du 



(24) 



where it can be shown that an outgoing wave (wave trans- 
porting angular momentum outward) is obtained when 
C = 0. This solution leads to the long "wakes" one ob- 
serves in compressible simulations of vortices. In particu- 
lar, in the limit x 3> kcg/qQ, we find that the wavefronts 
of the supersonic region are given by 



V = Vo 



2cs' 



(25) 



which produce a pattern similar to the wakes found in 
simulations (Fig.[5|. 

The connection of the supersonic and subsonic regions 
allows one to determine C and D knowing A and B. 
This can be done using asymptotic matching techniques 
based on Airy functions. Although the precise deriva- 
tion of these solutions is beyond the scope of this paper, 
we find that |C| ~ \D\ ~ max(|A|, \B\), i.e. the ampli- 
tude of the solution does not change significantly as one 



G 



crosses the sonic line. However, it should be noted that 
the wave amplitude in the supersonic region is signifi- 
cantly reduced compared to its amplitude at x = due 
to the exponential decay in the subsonic region. One can 
estimate this attenuation factor T with: 



Uy(±X s )\ 
\u y (x )\ 



'. exp 



= exp I 



A 2 H 



1/2 



du 



4qQkc s 



(26) 



This attenuation factor has a simple physical interpre- 
tation: it relates the vortex amplitude to the emitted 
density wave amplitude at a particular wavenumber. 
One can s how that the attenuation factor is minimum 
for k — q(2 — q)!!^ 1 where H is the disk thickness 

h = c s /n. 

It is clear that in the shearing-box the wave emitted 
in the inner side and the wave emitted in the outer side 
have the same amplitude (r + = r _ ). Because of this 
symmetry, one does not expect any migration. However, 
if one includes global effects such as curvature and the 
dependence of 51 as a function of radius, then an asym- 
metries might occur, resulting in vortex migration. 



background shear, is equal to the sound velocity. This 
corresponds to the large k limit in the local model. How- 
ever, this does not occur for small m. In fact (29) indi- 
cates that the distance of the sonic lines from the vortex 
increases as m decreases to the extent that when m = 1, 
the inner sonic line is out of the domain. This is because 
the phase velocity of the density waves increases as m 
decreases. WKB solutions very similar to the one pre- 
sented for the shearing box case can be derived and one 
can compute the attenuation factor as: 



r± = 



r.\u 



exp =F 



A 2 ( u ) 



1/2 



du 



(31) 

As it can be seen from this expression, the attenuation 
factor on the inner side (r _ ) is expected to be different 
from the attenuation factor on the outer side T + because 
A is not symmetric on both sides of corotation. One 
can easily check that |A 2 (ro + Sr)\ < |A 2 (ro — Sr)\ when 
Sr < rf — tq. Moreover, since |r+ — tq\ < \rj — ro|, we 
conclude that T + > r~, i.e. the outer wave always have 
an amplitude larger than the inner wave. 

This can also be seen by making the transformation 
from r to ip such that 



5.3. Cylindrical Geometry 

To extend our analysis to cylindrical geometry, we con- 
sider the case of a constant vortensity disk £ oc r~ 3 / 2 
with a constant aspect ratio c s cx r -1 / 2 . Furthermore, 
we assume the flow is barotropic and work in an inertial 
frame. We Fourier transform the velocity field in t and tp, 
defining u(t, tp, r) — J2ml dwexp[i(ujt — rn(p)]u(uj,rn,r). 
For a vortex corotating with the gas at r — ro we adopt 
lo = mQ(ro). We then perform the corresponding lin- 
earization of the equations to that described in the shear- 
ing box case. One eventually obtains: 



— —\lr—ru v +J\ (r)ru ip = 
ilr dr dr 



where 



A 2 (r) 



(T 2 (r) -n 2 



(27) 



(28) 



and er(r) — lj ~ mfl(r). This equation is very similar 
to the shearing-sheet problem. As in the shearing sheet 
case, we find two turning points rf on both sides of the 
corotation radius r defined as cr(r ) = 0. However, it 
should be noted that |r+ — r*o| < \r~ — ro|, i.e. the outer 
sonic line is closer to corotation than the inner sonic line. 



Indeed, imposing A( 



leads to: 



r t = r a I 



1 ± 



2/3 



(29) 



where e is the constant aspect ratio e 
m — > oo we have 



H/r. In the limit 



2/3 



ro 1± 



(30) 



Hence in this limit the sonic lines occur where the ve- 
locity relative to the vortex, that is associated with the 



^ = 1 
ro 



1 



cos ip, 



(32) 



so that the domain (rJ,fo) is mapped to (0,7r/2) and 
the domain (rJ,ro) is mapped to (tt/2,tt). The integrals 



inr d 



A 2 (u) 



1/2 



du 



(33) 



transform to 



lnT d 



-tt(-3±1)/4 
-7t(-1±1)/4 



2m 
"37 



(e 2 + ^)sin 2 ^ 



dip. 



(34) 

On account of the differing integration ranges for the 
outward and inward propagating waves the wave atten- 
uation factor is readily seen to be larger for the inward 
propagating wave. Although the calculation can be done 
in full analytically, it is adequate for our purposes to re- 
gard e and 1/m to be comparable small parameters and 
express the result correct up to order e. Then 



lnr d 



6 




(35) 



Thus the expected ratio of outward to inward wave flux 
arising from this asymmetry is 



exp 



— \ — me 
9 V m 



1 

me 



3/2' 



(36) 



As we will see later, this explains the inward migration 
found in simulations with a sim ilar profile. For example 
when e = 0.1 and m = 5, J36j) predicts an asymmetry 
comparable to that seen in Fig. [3] More generally, when 
m = 1/e we obtain a ratio ~ 1 — 2.5/m. 
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5.4. Vortensity gradient 

The above discussion on the local shearing box model 
assumed a constant vortensity and in this case the re- 
sponse to a symmetric vortex was found to be symmetric 
leading to equal inward and outward angular momentum 
flow rates and hence no net migration. However, if the 
constraint of constant vortensity is relaxed an asymme- 
try in the response occurs in the subsonic region of the 
flow even when the vortex is symmetric. To see this we 
consider a local shearing box model in which the back- 
ground density has the form E = const x cxp(— ax) with 
a being a constant. The vortensity of course varies as 
the inverse of E. 

As we are interested in the subsonic region we assume 
the linear response satisfies the anela stic condi tion V ■ 



(Eu) = 0. After eliminating E' from ([lT} and ([lSj) and 
then using the anelastic condition to eliminate u y , we 



obtain an equation for Q = v x y E in the form 



d 2 Q 

dx 2 



Q(h 



a 
'Ax 



(37) 



where k 2 



-a 2 / 4. We remark that for assumed unit 



radius of the center of the box, k y = m is the azimuthal 
mode number and a is the power law index occurring in 
E cx r~ a for a global disk. 



We are interested in solutions of (37) that decay expo- 
nentially for large |a;| and are matched as x — > 0. How- 
ever, we do dot require that derivatives match for x — > 0. 
This is because we assume that there is a nonlinear vor- 



tex core there that is not modeled by (37) which only 



describes the regions between the vortex core and the 
sonic lines. 

Without loss of generality we shall consider x > 0. 
The case x < can be recovered by letting a — > —a. 
An appropriate integral representation for the required 
solution follows fro m representations of 'Coulomb W ave 
Functions' given in 
may be written as 



Abramowitz & Stegun (19721 and 



C T(q)(2k)- q exp(-kx) 
1 

o \ 2fc 



x ■ 



t q exp(-t)dt. (38) 



where Co is an arbitrary constant and q = —a/ (6k). 
The above form may be used to connect the form of the 
solution near x = to that for large x giving 

Q(x = 0) = 1 -> Q 
~ T(q + l)(2kx)~ q exp(-kx) for x ->• oo (39) 

which may also be expressed as 

y(x = 0) = l^ 
y ~ exp(— kx — q\n(2kx) — "fq) for x — > oo. (40) 

Here, assuming that k is large, we have used T(q + 1) ~ 
— 7<7, where 7 is Euler's constant ~ 0.577. 

Recalling that q = — a /(6k), we find that the solution 
exceeds the pure exponential for large x > when a is 
positive. That means the outgoing wave is favored in 
this case. Conversely the ingoing wave is favored when 
a < indicating the potential for inward migration in 
this case. But, noting that kx ~ 1, this is found to 
lead only to a relatively weak modification of the wave 




Fig. 6. — Isocontour of vorticity for a vortex embedded in a 
vortensity gradient. The vortex is deformed by the stratification 
and gets flattened on the denser region (x < 0). 

propagation, except when \a\ is large and the vortensity 
gradient is very strong. 

However, it should be noted that the presence of a 
vortensity gradient leads to the modification of the vor- 
tex structure itself. In the above calculation, we have 
assumed the vortex core was symmetrical, i.e. u x (x — 
+e) = u x (x = — e) when e —> 0. However, in the pres- 
ence of a vortensity gradient, vortex solutions are not 
symmetric anymore. To show this effect, we have solved 
a set of anelastic equation in a shearing box including a 
radial vortensity gradient: 

V-(E (aj)t>) = ) (41) 

— +v-Vv = -VII - V(-qn 2 x 2 ) -2fixi), (42) 

In these equations, the vortensity gradient is imposed 
by the density profile Eo(x) and density/sound waves 
are filtered out. We have integrated numerically these 
equations assuming an exponential profile for the sur- 
face density E = Eoexp(— ax) using a spectral scheme. 
The initial conditions were chose n to be those appropri- 
ate to a Kida vortex (Kida||1981 ) of aspect ratio 8. This 



vortex then relaxes to a new quasi-stationary state corre- 
sponding to a vortex solution with a vortensity gradient. 
We show in Fig. [6] the resulting vorticity distribution in 
the case a — 4 (the length unit being the vortex ra- 
dial size). As it can be seen from this figure, the vortex 
core becomes asymmetric when a vortensity gradient is 
imposed. In particular, the perturbation extends for a 
greater distance on the less dense side and the vortex 
is "flattened" on the denser side. As a result, the less 
dense part of the vortex, being more extended, will ex- 
cite stronger density waves than the inner side, resulting 
in a new source of wave asymmetry. It should be noted 
that this effect is purely hydrodynamic and is due to the 
fact that the perturber (the vortex) "feels" and responds 
to the surrounding flow. Therefore, no corresponding 
effect exists in the case of a planet embedded in a disk. 

5.5. Comparison of vortex and planet migration 

We now summarize some important differences be- 
tween vortex migration and Type I planet migration 
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TABLE 2 

Comparing planet and vortex migration 



Effect 


Vortex 


Planet 


geometrical 


yes 


yes 


core asymmetry 


yes 


no 


corotation torque 


no 


yes 


pressure buffer 


no 


yes 



(see Table [2]). Although both are driven by asymmetric 
density waves, the excitation mechanism as well as the 
sources of the asymmetry are fundamentally different. 
First of all, it should be pointed out that although it has 
very similar properties, in contrast to the plan et m igra- 
tion problem, the asymmetry discussed in Sect. |5.3| does 
not involve any torque calculation. Instead, it appears 
as an attenuation in the linear propagation of density 
waves emitted by the vortex and it leads to a stronger 
outer wave. This is a major difference with the planet 
case in which asymmetric Lindblad torques acting at the 
sonic line are believed to drive planet migration. Fur- 
thermore, since we find the vortex always to move with 
the local gas velocity, there is no effect of the pressure 
buffer as discussed in the appendix. We find that wave 
asymmetries only weakly depend on the background den- 
sity profile, favoring the outer wave and therefore inward 
migration. 

However, unlike a planet, a vortex directly feels the 
disk and its associated density profile, which leads to an 
asymmetric vortex core as discussed in Sect. |5.4| in the 
case of a vortensity gradient. The resulting asymmetric 
wave emission can cither counteract or reinforce inward 
migration due to geometrical effects as discussed above. 

We also point out that since there is no gravitational 
interaction between disk and vortex, there are no horse- 
shoe orbits close t o corotation . In the planet case, 
the horseshoe drag (|Ward||199I I is an additional source 



of angular mo mentum exchange ( |Paardekooper fc Pa- 
paloizou|2009 ), that can be strong enough to counteract 
the wave to rque in non- isothermal disks ( Paardekooper 
et al. 2010 1. In the vortex case, we only have the wave 



torque. 

6. ANGULAR MOMENTUM TRANSPORT AND VORTEX 
MIGRATION 

It can be shown that the density waves d escribed above 
transport angular momentum outward (Papaloizou & 
Lin 1995). Therefore, if the waves emitted by a vor- 



tex are asymmetric then angular momentum has to be 
deposited or extracted in the vortex neighborhood. This 
in turn can be linked to vortex migration since moving a 
vortex from one location to another means some angular 
momentum was exchanged between these two locations. 
This argument shows that one can in principle relate the 
vortex migration rate to the wave emission asymmetry, 
without knowing precisely how the waves interact with 
the vortex structure. To quantify this effect more pre- 
cisely, let us start considering the angular momentum 
conservation equation: 



d t r{^(u v +rn{r)) + -d r r{F) 
r 







(43) 



where J- is the outward angular momentum flux J- = 
rT,u r (rfl(r) + u v ) = A/{2irr) (see equation Q) and (•} 
denotes a ip average. Let us assume a vortex of size 2s 



is located at r — ro and emits density waves asymmetri- 
cally. Integrating the above equation between the inner 
and outer radius tells us the vortex is losing angular mo- 
mentum if J-(ro + s) > J-(ro — s), i.e. if the outer wave is 
stronger than the inner wave, as in the constant vorten- 
sity disk described above. Next, to relate this effect to 
vortex migration, one has to compute the angular mo- 
mentum of a vortex. 

The angular momentum excess (compared to an un- 
perturbed disk) due to a single vortex is given by: 



J 



dr r 2 {8HrVL{r) + ^u^Ay (44) 



where S is the unperturbed density and the integration 
is computed between the inner and outer radius of the 
vortex. We denote the angular size of the vortex as A<p, 
which we expect to be ~ Hq/tq. It can be seen that the 
angular momentum of a vortex is made of two parts: a 
contribution due to the density fluctuation associated to 
the vortex, and a contribution due to the vortex rotation 
profile Utp. Assuming the vortex rotation rate is u>, the 
velocity can be approximated by u v ~ (r — tq)uj for ro — 
s < r < ro + s. The density fluctuation <5X can then 
be calculated assuming a geostrophic equilibrium in the 
vortex core. Using the radial component of the equation 
of motion we get: 



51(ro)u;s 2 



,5S 



(45) 



where we find as expected that anticyclonic vortices (u < 
0) are associated with high density regions. Combi ning 
the velocity profile and the density fluctuation in (44 1 
leads to: 



J ~ Eqws r Q Aip 



rgO(ro) 



1 



(46) 



where it can be seen that the density fluctuation con- 
tribution is larger than the rotation profile contribution 
by a factor (r/H) 2 3> 1. Therefore, a rather surprising 
result is that anticyclonic vortices are associated with an 
excess of angular momentum within the disk. If a disk 
region containing a vortex is losing angular m ome ntum 
by asymmetric wave emission, one deduce from ( 46 1 that 



the vortex can either migrate inward (reduce ro) or re- 
duce its size s. The procedure we have used to link wave 
emission to vortex migration does not allow us to directly 
compute the migration rate since the vortex size is also 
expected to change (in particular because s has to be 
kept smaller than H). However, it explains the inward 
migration observed in simulations of constant vortensity 
disks. Note that this asymmetry is cancelled out by the 
effects of a vortensity gradient in the sense that the mi- 
gration is expected to revers e for a large positive surface 
density gradient (see section 5.4). In practice migration 
reversal occurs for a <^ corresponding to a uniform 
surface density (see below and Fig. [7]). This is less ex- 
treme than for the case when the vortex is replaced by 
low mass protoplanet which would migrate inwards and 
require a <~ —1 when nonlinearit y is taken into account 
( Paardekooper & Papaloizou 2009 1 for migration reversal 
in isothermal disks. Even steeper profiles are needed in 
the linear migration regime. This because the pressure 
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Fig. 7. — Migration speed in disks with Ho = O.lro and the 
same initial perturbation as in Fig. [l]for different surface density 
profiles. 

buffer associated with the fact that an orbiting planet 
moves at different angular velocity to that of the disk 
gas when there is a pressure gradient, which would cause 
a greater asymmetry in favor of the outer wave, in fact 
does not operate for a vortex (see the appendix). Since 
the vortex is tied to the disk, it will always move with 
the local gas velocity. 

7. GLOBAL COMPRESSIBLE SIMULATIONS 
7.1. Migration for different surface density profiles 

We have measured the migration speed in isothermal 
disks with Hq = O.lro for various background density 
gradients; the results are depicted in Fig. [7J The mi- 
gration speed is approximately proportional to —a, with 
outward migration happening for surface density profiles 
that increase outward. As indicated above a constant 
surface density disk shows almost no vortex migration. 
For such a disk, the geometrical effects favoring inward 
migration (see Sect, pi) are almost exactly balanced by 
the asymmetry of the vortex, which favors outward mi- 
gration. The vortex asymmetries make vortex migration 
more strongly dependent on a than is the case fo r lin- 
ear Type I planetary migration ( Tanaka et al.|2002[ ) that 
requires a <^ —2.5 for migration reversal. 

The asymmetry in the vortex can be seen for exam- 
ple in the Fourier components of the radial velocity. In 
Fig. [8j we show the m = 10 component for a constant 
vortensity disk (black curves) and a constant surface den- 
sity disk (grey lines). The approximate size of the vor- 
tex is indicated by the vertical dotted line. The con- 
stant vortensity disk gives rise to a vortex that is almost 
symmetric for small as expected. The constant sur- 
face density disk shows a stronger response in the inner 
disk (dashed curve) , which favors outward migration and 
works against the geometrical effects of Sect. [5] From 
Fig. [7] we see that for this disk, these competing effects 
almost cancel each other, resulting in almost no migra- 
tion for a = 0. We comment that for a disk with a 
temperature gradient (through a spatially varying, but 
constant in time, sound speed), this neutral state corre- 
sponds to a disk with constant pressure. 

7.2. Migration and disk thickness 
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Fig. 8. — Variation of the m = 10 component of the radial ve- 
locity with distance to the vortex x = (r — r VO rt ) /r V ort , in an 
isothermal disk with Ho = O.lro. The regions inside and outside 
of the vortex are represented by dashed and solid curves, respec- 
tively. Two density profiles are considered: a = 3/2 (black curves) 
and o = (grey curves). The vertical dotted line indicates the 
approximate extent of the vortex. 




Fig. 9. — Top panel: wave action due to a vortex of size Ho/2, 
and initial velocity perturbation 3c 3 /4, in disks with a = 3/2 for 
different Ho- Circles indicate the wave action inside of the orbit 
of the vortex, triangles the wave action outside the orbit of the 
vortex, and diamonds the difference between the two. Curves are 
fits through the data <x H^ (for Ai and A ), and oc Hq (for AA). 
Bottom panel: corresponding migration speeds. Squares indicate 
numerical results, the solid curve if a fit through the data <x Hq. 

Until now, we have focused on a disk with Hq = O.lro. 
Since we need to be able to resolve the vortex, which 
is of size H /2 typically, it is numerically convenient to 
work with relatively thick disks. However, we expect 
the angular momentum transport, and therefore the mi- 
gration speed of the vortex, to depend on Hq. From 
equation |9]), we expect vortices of size ~ Hq and veloc- 
ity perturbations ~ c s , that A oc Hq. Since the most 
important contributions to the angular momen tum flux 
come from m ~ ro/i?o> we see from equation (36) that 
the ratio of wave actions outside the orbital radius of 
the vortex and inside the orbital radius of the vortex 
is A Q /Ai w 1 — CHo/ro, with C a numerical constant 
of order unity. This means that the action difference 
A A = Ai — A Q ex Hq. This equals the rate of angular 
momentum transfer between vortex and disk, and there- 
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Fig. 10. — Top panel: wave action due to vortices of different 
sizes, with initial velocity perturbation 3c s /4, in disks with a = 3/2 
and Hq = O.lro. Circles indicate the wave action inside of the 
orbit of the vortex, triangles the wave action outside the orbit of 
the vortex, and diamonds the difference between the two. Bottom 
panel: corresponding migration speeds. 

fore governs vortex migration. 

In the top panel of Fig. [9] we show the measured wave 
action inside the orbit of the vortex (circles) , outside the 
orbit of the vortex (triangles), and their difference for 
different values of Hq. The size of the initial perturba- 
tion is a fixed fraction (1/2) of Hq, and the magnitude of 
the velocity perturbation is a fixed fraction of the sound 
speed (3/4). With these initial conditions, we expect 
above scalings to hold approximately. We have overplot- 
ted fits through the numerical results oc Hq (for Ai and 
A a ) and oc Hq (for AA), showing that indeed these scal- 
ings hold for a wide range of Hq , even in cases where the 
disk can no longer be called 'thin'. 

The angular momen tum of the vortex scales as Hq, 
according to equation ( 46 1 . Combined with an angular 



momentum transfer rate oc Hq as derived above, we de- 
rive a migration time scale oc Hq 2 , or a radial velocity of 
the vortex oc Hq. This is the relation we explore in the 
bottom panel of Fig. [9j where we show the magnitude of 
the migration velocity of the vortex for different values of 
Hq. The solid curve is a fit through the numerical data 
oc Hq, showing that again the scaling derived shows good 
agreement with the numerical simulations. Note that un- 
like for the case of planet migration, vortex migration is 
faster in thicker disks. This is because the size of the 
vortex scale with Hq, so that for all Hq we get a similar 
perturbation. This would amount to considering planets 
of higher mass in thicker disks. 

Even though vortex migration slows down for thinner 
disks, it is still very substantial at H = 0.05r , a value 
that is though to be appropriate for protoplanetary disks. 
From the measured radial velocity, we derive a migration 
time scale of ro/\v r \ = 1700 orbits. This is very short, 
but we point out that these results were obtained for a 
strong vortex with a size comparable to Hq. Below, we 
discuss how migration depends on the size of the vortex. 

7.3. Effects of vortex size and strength 

Vortex migration is driven by perturbations at the 
sonic line, located roughly at a distance 2H /3 away 
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Fig. 11. — Top panel: wave action due to vortices of different 
aspect ratio x, with radial size Ho/2, in disks with a = 3/2 and 
Hq = O.lro. Circles indicate the wave action inside of the orbit 
of the vortex, triangles the wave action outside the orbit of the 
vortex, and diamonds the difference between the two. Bottom 
panel: corresponding migration speeds. 

from the center of the vortex. Away from the core, per- 
turbations due to the vortex de cay exponentially in the 
subsonic region (see equation (26)). We therefore ex- 
pect smaller vortices to show less migration. Note that 
a smaller radial vortex size s can be compensated for in 
principle, as far as migration is concerned, by a stronger 
initial velocity perturbation. Unfortunately, the regime 
s -C Hq is difficult to study numerically because of com- 
putational constraints. We therefore restrict ourselves to 
s > Hq/10 to study the effect of vortex size on migration 
We show numerical results for Hq 
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O.lro in Fig 

Note that the radial size s quoted in the figure is actually 
the size of the initial perturbation rather than the actual 
vortex size. We have always found the two to agree well. 
We see that indeed smaller vortices migrate slower, with 
a factor 50 between a vortex of size Hq/10 and a vortex of 
size Hq/2. For s > Hq/2, migration is roughly constant, 
but for the larger cases the vortex extends beyond the 
sonic line, a situation that is difficult to maintain. 

It is difficult to provide simple scalings in this case, 
since for most of our results the vortex core occupies a 
substantial fraction of Ho. We have checked that similar 
results are obtained for different values of Hq. Therefore, 
even a vortex of size Hq/10 with Hq = 0.05ro would have 
a migration time scale of ~ 30000 orbits only. 

As mentioned above, a stronger velocity perturbation 
will induce stronger wave emission, and therefore faster 
migration. For a given radial vortex size, a stronger ve- 
locity perturbation will lead to a vortex that is more 
circular, while weaker vortices have a larger aspect ratio 
X- To study the influence of the strength of the vortex, 
we have varied the initial velocity perturbation at a fixed 
radial size of the vortex. The aspect ratio of the resulting 
vortex was measured by considering contours of constant 
vortensity. 

The resulting wave action and migration speeds are 
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Fig. 12. — Top panel: surface density profile used to study the 
stalling of vortices. Bottom panel: the corresponding profile of 
a = — dlog S/dlog r. 
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Fig. 13. — Top panel: evolution of the strength of the vortex, 
defined by the minimum in the relative vortensity perturbation. 
Middle panel: the local value of a. Bottom panel: orbital evolution 
of the vortex. 

displayed in Fig. [TTJ As expected, stronger vortices 
migrate faster. This is not only due to a larger per- 
turbation amplitude in the velocity, but also due to a 
change in the perturbation spectrum. Since the attenu- 
ation factor peaks at to ~ to/Hq, the perturbation am- 
plitudes around this value of to play a major role in de- 
termining the migration speed. A vortex of aspect ratio 
unity and radial size H gives rise to large amplitudes 
at to ~ ro/Ho, and therefore to strong migration. In- 
creasing the aspect ratio shifts the peak of the vortex 
spectrum to lower values of to, which are less effective 
in inducing angu lar momentum transfer. This is what is 
observed in Fig. [TTJ Migration slows down by an order 
of magnitude when the vortex weakens from \ = 2 to 
X = 3.5. 

7.4. Vortex trapping at surface density maxima 

The strong dependence of migration rates on the back- 
ground surface density profile opens up the possibility 
of halting vortex migration at special locations in the 
disk, where the local gradient of surface density is small. 
This is explored using a density profile as depicted in 
the top panel of Fig. [12] The corresponding profile of 
a = — d log £ jd log r is shown in the bottom panel of Fig. 
12 From the discussion above we expect a vortex to mi- 



grate inwards outside r = rn, where a > 0, and outward 
inside r — vq, where a < 0. In this configuration, r = 7'o 
is an equilibrium radius for migrating vortices. 

We have put in a strong vortex, with initial veloc- 
ity perturbation amplitude of 3c s /4 over a scale Hq, at 
r/r = 1.5 in a disk with H Q — 0.1r . With such a strong 
vortex, we make sure it does not dissipate before reaching 



13} we show the time evo- 
le vortex (bottom panel), 



the equilibrium radius. In Fig. 
lution of the orbital radius of tf 

the strength of the vortex S v (defined by the minimum in 
the relative vortensity perturbation, top panel), and the 
local value of a (middle panel). From the top panel we 
see that the vortex gets weaker with respect to the back- 
ground flow, partly because of numerical diffusion, partly 
because of the migration itself. Since the vortex, being 
a negative perturbation in vortensity, is moving into a 
region of lower vortensity, the relative perturbation is 
expected to go down, as observed in Fig. [13] Results 
for two different resolutions are shown in Fig. [13] the 
evolution of S v that appears similar in both cases is due 
to migration, while differences are most likely caused by 
numerical diffusion. 

From the bottom panel of Fig. [13] we see that mi- 
gration starts off fast, slows down over time, and almost 
stalls at r — ro- This slowing down of migration is com- 
pletely due to the decrease in a as depicted in the middle 
panel of Fig. 13 At r/r — 1.5, a « 4, giving rise to 
fast inward migration. As the vortex approaches the sur- 
face density maximum, a approaches zero and migration 
slows down. Since we obtain essentially the same result 
for both resolutions, we conclude that the halting of in- 
ward migration is due to r — r$ being an equilibrium 
radius rather than due to the vortex dissipating. This 
therefore provides an example of vortex trapping at a 
surface density bump. 

8. DISCUSSION AND CONCLUSIONS 

We have identified a mechanism for angular momen- 
tum exchange between a vortex and the surrounding 
gaseous disk, which leads to orbital migration of the 
vortex. The vortex excites spiral density waves inside 
and outside its orbit that will in general be of unequal 
strength. Geometrical effects favor the outer wave, while 
a background vortensity gradient in the disk leads to an 
asymmetric vortex, favoring the wave where the vorten- 
sity is highest. For a constant surface density disk, the 
two effects nearly cancel. 

Whenever the waves are asymmetric, there is an asso- 
ciate angular momentum change in the region around the 
vortex. The vortex can either shrink, or change its or- 
bital radius. Numerical simulations indicate that vortex 
migration can be very fast, especially for strong vortices 
embedded in thick disks, for which we found a migration 
time scale of a few 100 orbits for H = 0.1r . Migration 
slows down rapidly for vortices with size s -C Hq. Im- 
portant differences between vortex migrati on a nd Type I 
planet migration are summarised in Sect. 5.5 and Table 



m 

A few simplifications were made in the analysis. First 
of all, we have not solved the energy equation, working 
with a prescribed fixed sound speed. We do not expect 
things to change when including the full energy equation. 
Preliminary simulations show rapid vortex migration in 
this case as well, with a constant pressure disk again 
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being the neutral state. Migration rates are comparable 
to the isothermal case. 

We have neglected the self-gravity of the disk. Since 
the vortex is associated with a surface density maximum, 
it will act as an embedded mass as well as an obsta- 
cle. The gravitational interaction of vortex and disk can, 
completely analogue to the case of an embedded planet, 
lead to additional migration. However, unless the disk is 
extremely massive, this effect will be small. This is be- 
cause a vortex of roughly size Hq perturbs the flow like a 
planet for which the Hill sphere is comparable in size to 
Hq, or typically a Saturn-mass planet for H — 0.05r . 
On the other hand, for a Minimum Mass Solar Nebula 
disk, the mass contained in the vortex is less than an 
Earth mass. This is one of the reasons why vortex migra- 
tion can be much more efficient than planet migration: 
we have an Earth-mass object, emitting density waves 
like Saturn. 

We have worked in the two-dimensional approxima- 
tion. Although vortices ar e in general unstable in 3D 



(Lesur & Papaloizou 2009), they can survive if an ex- 



citati on mechanism operates (e.g. |Lesur fc Papaloizou 
20101. The two-dimensional approximation is then valid 



as long as the vertical size of the vortices produced is 
comparable to Hq, since one can then work with verti- 
cally averaged quantities. Note however th at vortices can 



also i nduce large scale 3D circulations (e.g. Meheut et al 



20101, in which case the two-dimensional approximation 
might not be totally valid. 

We have considered a single vortex only. Although a 
detailed description of the interaction of multiple vortices 
is beyond the scope of this work, we note that in a 2D 
disk, vortices quickly merge, forming a larger vortex that 
subsequently migrates faster than the original vortices. 

We have not considered the effects of turbulence in the 
disk. Our results would therefore apply to a region of the 
disk where turb ulence is almost absent, like a dead zone 
(Gammie 1996). A thorough discussion of the possible 
effects of turbulence on vortex migration is beyond the 
scope of this paper, but we note that if a vortex can sur- 
vive in a turbulent disk, it would be subject to the same 
wave torque as discussed in this paper. However, since 
the wakes are now generated in a stochastic flow, there 
will be a stochastic component in the torque. For the 
case of an embedded planet, where the interaction with 
the turbulenc e is gravitational, this is kno wn as stochas- 
tic migration ( |Nelson fc Papaloi zou 2004]). 

The efficiency ot vortex migration raises some impor- 
tant questions about the role of vortices in planet for- 
mation. It is well-known that anticyclonic vortices can 
trap dust particles, and may therefore be efficient sites 
of planet formation. Now unless vortices are formed at a 
special location in the disk that is neutral to vortex mi- 
gration (i.e. a region of constant pressure), the collected 
solids will rapidly migrate inward with the vortex, all the 



way to the inner edge of the disk if there is no pressure 
maximum in between. Therefore, if vortices play a role 
in planet formation through accumulating dust particles, 
all these solids will end up either at a special location in 
the disk where vortex migration stalls, or at the inner 
edge of the disk. 

In the absence of any stalling location, the life time of 
a vortex is basically given by its migration time scale, 
which is only a few thousand orbits for a vortex of size 
Ho/2, aspect ratio \ — 2 and Hq — 0.05ro. This pro- 
vides a strong time scale constraint for possibly forming 
long-period (i.e. orbiting at several AU) planets inside 
a vortex, since they need to decouple from the vortex 
before reaching the inner edge of the disk. 

Vortices produced by a Rossby wave instability 
( |Lovelace et al||1999[ ) are formed near pressure maxima, 
and therefore find themselves in a migration neutral lo- 
cation from the start. However, this instability requires 
strong variations in surface density within one pressure 
scale height, and it is not clear if such a configuration can 
be formed naturally in the first p lace. One possibi lity is 
the outer edge of the dead zone (Lyra et al. 20091), but 



this would require a very sharp transition in resistivity. 

Vortices appear to be naturally produced in g lobal sim- 
ulations of the MRI ( Fromang & Nelson 2005 1 , and can 
also grow in regions of a negati ve entropy gradient by the 
SBI ( |Lesur fc Papaloizou|2010 ). In both cases, these vor- 
tices wolIH - b^~s1lD3ecTToiriigration. This way, the outer 
edge of a dead zone can harbour vortices even in the 
absence of strong gradients locally, since even a smooth 
pressure maximum will trap incoming vortices that were 
generated in the outer disk by either the MRI or the SBI. 
Note that vortex migration makes it more difficult to sus- 
tain the SBI locally. Since the SBI is subcritical, a finite 
amplitude perturbation is required to start it. Once vor- 
tices of size ~ Hq form, they quickly disappear into the 
inner disk, after which another perturbation is needed to 
get a second generation of vortices. 

In conclusion, vortex migration may be a very effective 
way of collecting solids in special locations in the disk, 
such as pressure bumps. Plan et formation can the n pro- 
ceed rapidly by coagulation (Brauer et al. 2008), even 
if the original vortices do not survive. If vortices form 
readily in protoplanetary disks, this suggests that pres- 
sure maxima will be preferred sites for planet formation. 
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APPENDIX 
PRESSURE BUFFER 



We have found the vortex to always orbit at the local gas velocity. This is a different situation to that appropriate 
to a body that is decoupled from the disk flow, such as a planet, on which pressure forces do not act, and therefore 
would move with the local Keplerian angular velocity whereas the gas moves with a sub-Keplerian velocity. For the 
constant vortensity constant aspect ratio disk model discussed in[5j such a body located at r = tq would rotate with 
the Keplerian angular velocity \Q(r ), where A = 1 + 5e 2 /4. This velocity difference leads to a significantly increased 
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asymmetry in the wave amplitudes that is said to be a consequence of the pressure buffer that slows the rotation of 
the gas relative to the orbiting object for the disk model considered here. 

It is a simple matter to repeat the calculation in Sect. [5] for the wave attenuation factors appropriate to this case 
after making the replacement u = m\Q(ro). In this case we make the transformation 



= 1 



'0 



A - 1 + e 2 + 



COS 1p, 



(Al) 



so that the integrals determining them become 



2m 



(A-l + e 2 + ^)sin 2 ^ 



3e 



V\ + yj \ - 1 + e 2 + ^ COSTp 



di>. 



The ratio of outward to inward wave flux corresponding to ( 36 1 is then 

2 

= exp 



9\l m 



9 



-me H 

4 me 



3/2' 



(A2) 



(A3) 



When m = 1/e this leads to a ratio 1 — 5.2/m, corresponding to an asymmetry about twice as large as given by (36) 



Therefore, the pressure buffer is an import ant addition al source of asymmetry in the case of an embedded object on 
a Keplerian orbit such as a planet (see also Ward|[l997 1 . 
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